Classical invariants and the quantization of chaotic systems 
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Long periodic orbits constitute a serious drawback in Gutzwiller's theory of chaotic systems, and 
then it would be desirable that other classical invariants, not suffering from the same problem, 
could be used in the quantization of such systems. In this respect, we demonstrate how a suitable 
dynamical analysis of chaotic quantum spectra unveils the fundamental role played by classical 
invariant areas related to the stable and unstable manifolds of short periodic orbits. 
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The correspondence between quantum and classical 
mechanics has been a topic of much interest since the 
beginning of the quantum theory, and more recently in 
relation to quantum chaos 0, Q . The question involves 
elucidating the classical objects and properties on which 
to impose quantum restrictions, this being at the heart 
of every semiclassical theory. 

Very early, Ehrenfest noticed Q the importance of 
classical adiabatic invariants, such as the action, in the 
quantization of dynamical systems. Later, Einstein Q 
realized that the proper arena to perform this quantiza- 
tion for integrable motions are invariant tori He also 
remarked that this theory is not applicable to chaotic 
motions, due to the lack of supporting invariant classi- 
cal structures. After that, dynamical invariants are re- 
garded as the geometrical objects upon which reason- 
able semiclassical theories of quantum states should be 
constructed. Concerning the associated properties to be 
used, those which are canonically invariant seem to be 
the natural choice, since they render descriptions inde- 
pendent on the coordinate system. 

Keeping within this scheme, Gutzwiller took in the 
1970's a new route, and chose periodic orbits (POs), and 
their individual properties (actions, Maslov indices and 
stability matrices), leaving aside others (corresponding 
manifolds), as the quantizable invariants In this way, 
he developed a semiclassical version of the quantum me- 
chanical Green function, that has become the cornerstone 
of the semiclassical quantization of chaotic systems. The 
resulting density of states appears as the sum of contri- 
butions of all POs of the system, and their repetitions. 
The phase of each contribution is the action (symplectic 
area) along the orbit (divided by %), and the amplitude 
is proportional to its stability. Unfortunately, this the- 
ory suffers from a serious computational problem in the 
exponential growth of the number of contributing orbits 
with the Heisenberg time, £h = 2nTip(E), with p{E) the 
energy density. This have precluded its use except in very 
special situations Q. In this respect, it is worth empha- 



sizing that Gutzwiller's summation formula can be used 
in two opposite ways. In the direct route, it can be fed 
with classical information to predict quantum eigenval- 
ues. Or alternatively, it can be used in an inverse way to 
extract classical magnitudes from the eigenvalues spec- 
trum |8( . Curiously enough, these two operations are not 
equivalent, in the following sense. When applied in the 
direct way, one needs to included longer and longer POs 
(with periods of the order of in) to predict higher energy 
eigenvalues. However, when applied backwards, for ex- 
ample, by Fourier transforming the eigenvalues spectrum 
of a chaotic billiard, the periods (or other properties) of 
only short POs, with values up to the magnitude of the 
Ehrenfest time, £e 0- are obtained [Hj- This asymmetry 
is not fully understood yet, and raises fundamental ques- 
tions about our present understanding of the quantum 
mechanics of chaotic systems: are long POs (with peri- 
ods of the order of in) relevant?, or its inclusion in the 
theory of Gutzwiller is only a drawback?, finally, respon- 
sible for the unreasonable computational effort involved 
in the semiclassical computation of physical magnitudes. 

In this Letter, we address this issue, by investigating 
the inverse route in a non-standard way. By explicitly 
including the dynamics of short POs (the only relevant 
in this route) in the Fourier transform process, we de- 
velop a method, relying only on purely quantum infor- 
mation, able to extract the pertinent associated infor- 
mation from the actual full quantum dynamics of very 
chaotic systems. We have found conclusive evidence that 
the corresponding quantum spectrum contains informa- 
tion about collective invariant objects associated to short 
POs, namely, the homoclinic and heteroclinic areas en- 
closed by their stable and unstable manifolds. This im- 
plies some sort of interaction between periodic structures, 
that can play a role equivalent to that of long POs in the 
Gutzwiller formula. 

To gauge the dynamical interaction between two POs, 
A and B, in a quantum sense, we propose the use of the 
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cross correlation function 
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where U(t) is the time propagation operator, and </>a,b 
are suitable functions associated to the POs, whose na- 
ture will be discussed later. In our formula, the sec- 
ond part of the bracket follows the evolution of the non- 
stationary function associated to one of the POs, and 
the application of the bra extracts the information rela- 
tive to the other PO contained in it, thus filtering out (at 
least to some extend) the rest. By choosing a correlation 
function as our indicator, we have the same information 
as in the corresponding spectra, but recast with a more 
dynamical perspective. 

A natural choice for <^a and (f>B are wavefunctions liv- 
ing in the vicinity of the corresponding POs. These 
functions are constructed very efficiently, either by dy- 
namically averaging over the short time dynamics of the 
associated PO 11], or by minimizing the energy disper- 
sion in a suitable basis of transversally excited resonances 
|l2|. In this work, we use scar functions as defined in 
Ref. 0. These functions are highly localized in energy 
around some mean values satisfying a Bohr-Sommerfeld 
type quantization rule 
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where S(E) is the dynamical action at energy E, v the 
Maslov index, and n an integer number. 

In order to study the previous ideas, we choose a parti- 
cle of mass 1/2 enclosed in a fully chaotic desymmetrized 
stadium billiard of radius r — 1 and area 1 + 7r/4, with 
Dirichlet conditions on the stadium boundary and Neu- 
mann conditions on the horizontal and vertical symme- 
try axis [see the inset in Fig.n(a)]. For this system, the 
action takes a simple semiclassical relation in terms of 
the mean wavenumber, k, and the length, L, of the PO; 
namely, S(E)/H — kL. 

In our numerical study, we consider the horizontal 
(A) and V-shaped (B) POs with lengths La = 4 and 
Lb = 2(1 + V2), respectively [see the inset in Fig.Q(a)]. 
Let 4>a be the scar function for A with mean wave num- 
ber kA [obtained from Eq. J5J], and 4>b the correspond- 
ing for B, with the mean wave number &b closest to &a- 
We focus in the cross correlation function as defined in 
Eq. and accordingly, we present in Fig. Q (a) C(t) 
for kA = 155.116. The most striking feature in the plot 
is the totally different behavior exhibited by the correla- 
tion, below and above times of the order of the Ehrenfcst 
time, t~E, (the actual value of which has been marked 
with an arrow in the figure). For short times, the corre- 
lation function increases monotonically from zero up to 
the first maximum. This maximum appears at approxi- 
mately twice the value of tE, point at which interference 
starts to be relevant. After that, other maxima appear, 
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FIG. 1: (a) Cross correlation function between the horizontal 
(A) and the V-shaped (B) periodic orbits shown in the inset, 
for a value of the wavenumber of &a = 155.116. The time t is 
measure in units of the Heisenberg time, 
(b) Maximum value of the cross correlation function in the 
interval [0, to] as a function of the wavenumber, for to = t E /4] 
(lower curve) and to = £h (upper curve). The mean decreas- 
ing tendency is indicated in dotted line. 



and the behavior of the correlation gets very complex for 
times of the order of the Heisenberg time, in, which is 
equal to one in the units system used by us. 

To further characterize the interaction between POs, 
some representative dynamically meaningful magnitude 
along the spectra should be defined. For this purpose, 
we take the maximum of C(t) in the time interval [0,to] 



C max (t o , kA ) 



c{C(i); for < * < t }, 



where the dependence on kA has been explicitly included. 
[For instance, in the case of to = tn, this maximum ap- 
pears marked with an asterisk in Fig. [I] (a)]. Obviously, 
other quantities, such as for example the integral of C(t) 
in the interval, can be used as the representative magni- 
tude; our conclusions, however, are independent of this 
choice. In Fig. (b) we show C max (to,fcA), as a func- 
tion of kA, for two values of the maximum time to, that 
have been taken equal to £e/4 (lower curve) and in (up- 
per curve). As can be seen, both functions decay as A; a 
increases, while oscillating at the same time with a dom- 
inant frequency. To analyze the frequency dependence of 
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FIG. 2: (a) Fourier transform of the maximum cross correla- 
tion functions in Fig. 0(b): to = te/4 (full line), and to = tu 
(dashed line). 

(b) Rescaled intensity of the Fourier transform for to = tn in 
part (a) after the big peak has been removed. 

these functions, we have Fourier transformed them, after 
the signal has been properly prepared by eliminating the 
decaying tendency (dotted line). The resulting spectra 
are shown in Fig. [21 (a) in full and dashed line, respec- 
tively. As can be seen, they both appears dominated by 
only one peak, at values of the action: S = 0.227 for 
to = *e/4, and S = 0.828 for t = t R . Notice that S has 
units of length due to the fact that the total linear mo- 
mentum of the particle has been set equal to one. From 
the discussion above, our aim is to correlate these two 
peaks with invariant classical structures related to the 
chosen POs. Taking into account the numerical values of 
their positions, the first peak (labelled HE in the figure) 
can be assigned to the heteroclinic area, Sab, enclosed 
by the stable and unstable manifolds emanating from the 
fixed points associated to POs A and B. This region is 
shown, shaded with horizontal lines, in the phase space 
portrait of Fig. |3 which illustrates the classical struc- 
tures relevant to our work. When calculated, the corre- 
sponding area is 0.22540, value that agrees extremely well 
with that numerically found for Sue- Moreover, accord- 
ing to the results reported in Ref. [l3| for a more abstract 
case, this heteroclinic area is related to the semiclassi- 
cal Hamiltonian matrix element between scar functions 
through the relation 

|(0b|^|0a)| 2 <xcos(Sab k), (3) 

which presents the same oscillating frequency found by 
us. This fact can be taken as a further confirmation of 
our previous assignment. 

Using the same kind of arguments as before, the peak 
on the dashed curve of Fig. [21 (a), corresponding to to = 
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FIG. 3: Phase space portrait in Birkhoff coordinates of the 
classical structures relevant to our calculations. The fixed 
points marked with (A) and (B) correspond to the labeled 
turning points of periodic orbits A and B shown in the inset to 
Fig-0(a) . The unstable and stable manifolds of the horizontal 
PO are represented in full and dashed line, respectively, and 
those for the "V" -shaped PO in dotted line. The shaded 
regions corresponds to the heteroclinic (horizontal lines) and 
homoclinic (vertical lines) referenced in the text. 

in, can be assigned to the difference in length between 
POs A and B. This quantity amounts to 0.828426 . . ., 
again in excellent agreement with the value found nu- 
merically. The existence of this peak, which is associated 
to the difference in actions between orbits A and B [see 
Eq. @] , reflects the strong dependence of the interaction 
with &a — &b (which can also be related with the energy 
separation between the resonant levels corresponding to 
<Aa.,b)- This effect has been further confirmed by ana- 
lyzing the oscillatory behaviour of this wavenumber dif- 
ference, which turns out to be the same exhibited by 
Cmax(tH,^A), although these two functions appear with 
opposite phases. 

At this point it is worth to emphasize that the differ- 
ent origin of the two peaks reflects the existence of two 
regimes in the cross correlation function (JJ, with the 
corresponding transition taking place at t ~ t-&. These 
two regimes can be easily understood in terms of the 
Fermi golden rule, since C max (ici) &a) is in fact a measure 
of the probability transition between the resonant states 
4>a and 4>b ■ Accordingly, the behavior of C max is given by 
the competition between two factors: the square of the 
coupling matrix element and the separation in energy of 
the corresponding levels. 

Let us now consider other components in the spectra of 
Fig. [21 (a). To observe them more clearly, we calculate the 
rescaled intensity that is obtained after the biggest peak 
has been removed. When this is done for the two plot- 
ted curves only the results corresponding to the dashed 
one are stable against local variations of t$. They are 
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shown in the lower part of the figure. As can be seen, 
many different contributions appear, all with a compara- 
ble order of magnitude. Among them we have focussed 
in the highest one (labelled HO), as the most interest- 
ing. When the value of Sho is compared to the relevant 
classical invariants of our problem (see Fig. we find 
that it matches remarkably well with the homoclinic area 
enclosed by the stable and unstable manifolds emanating 
from A; this region appears shadowed with vertical lines 
in Fig. |2| The assignment is also supported by the work 
by Ozorio de Almeida |14| , which has shown how homo- 
clinic motions can be quantized, using the corresponding 
invariant manifolds. We believe that all (or most) re- 
maining peaks in the spectrum presented in the lower 
part of Fig. can be interpreted in the same way, using 
different homoclinic and heteroclinic regions correspond- 
ing to the same POs. Actually, we have succeded in as- 
signing the two peaks located to the left of that at Sho > 
however, a full description of the procedure is deferred to 
a future publication. 

In summary, some aspects concerning the role of short 
POs in Gutzwiller's summation formula, the cornerstone 
for the semiclassical quantization of chaotic systems, have 
been analyzed. For this purpose, only purely quantum 
information has been used, in order to obtain relevant 
classical invariants for a semiclassical theory of quantum 
chaos. We have found evidence that the short POs, and 
their associated heteroclinic and homoclinic intersecting 
areas, are relevant contributions to the spectra. Further- 
more, our results provide a semiclassical interpretation on 
how structures localized along POs interact pjj- This in- 
teraction is found to be given as the competition of two ef- 
fects: the energy separation between scar functions, and 
the squared coupling matrix clement between them, fol- 
lowing the celebrated Fermi golden rule. Also, we have 
shown how the values of the corresponding parameters 
can be theoretically evaluated a priori. Finally, our re- 
sults point out to the possibility of constructing compu- 
tationally tractable alternatives to Gutzwiller's theory, 
in which the long POs arc substituted by the interaction 
between short POs. In this respect, the steps given in 
Ref. ^3 provide a suitable frame in this direction that, 
together with a deeper understanding of the interaction 
between POs, can be a bridge towards a fully satisfactory 
semiclassical theory of chaotic systems based on short 



POs. This theory would be interesting not only from a 
fundamental point of view, but also for its application, 
for example, to nanotechnology Q, where it has been re- 
cently shown useful, for example, in the study of electron 
transport in mesoscopic devices |l7|. 
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